globals;
[priors,mmus,rr] = gmmparams(10,5);
[X,labels] = gmmrnd(priors,mmus,rr,1000);
f=figure;
gmmplot(X,mmus,rr,f);

